32nd International Cosmic Ray Conference, Beijing 201 1 



o 

> 
O 

:^ 

On 
W 

X 

6 



> 

in 



X 
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Abstract: A few Galactic objects are known to be variable sources of photons with energies above 100 GeV. These 
systems are mostly binaries, where variabihty can often be connected to the orbital period. Particle acceleration and 
gamma-ray production processes in binaries are generally not well understood. We present here an overview of VERITAS 
observations of binary systems at very-high energies (VHE) with emphasis on LS I -1-61 303 and 1A0535-I-262. The results 
are discussed in the context of recent observations of RXTE, Swift and Fermi LAT at X- and Gamma-ray energies. 

Keywords: acceleration of particles ; binaries; gamma rays: observations - individual (lA 0535-1-262; LS I -1-61 303) 



1 Introduction 

Several high-mass X-ray binaries have been detected in the 
VHE energy range in the past years. These binaries con- 
stitute the only known class of galactic objects with vari- 
able point-like VHE emission. The class currently contains 
three members: PSR B1959-63/LS 2883 [6], LS 5039 [7] 
and LS I +61 303 [9, 2]. A fourth TeV binary might be the 
unidentified VHE source HESS J0632+057 [8], a variable 
gamma-ray source located in the Monoceros region. 

The variable emission from VHE binaries is likely con- 
nected to changes in physical parameters associated with 
the orbital movement. This makes each of these systems 
unique, since orientation, sizes, and eccentricity of the or- 
bits vary widely. Orbital periods range from a a few days in 
LS 5039 to several years in PSR B 1959-63. The changes 
during an orbit are not of geometrical nature only, beside 
distances and orientation of the objects towards each other 
and the observer are variations in photon and matter den- 
sities around the acceleration sites the main causes for the 
variable emission. 

PSR 1259-63 is currently the only VHE binary with a 
clearly identified compact object - a 48 ms pulsar orbiting 
a Be star in an eccentric orbit. No pulsed emission at any 
wavelength has been detected from the binaries LS 5039 
and LS I + 61 303, although pulsed signals might get ab- 
sorbed in the dense photon field of the nearby bright com- 
panion stars (see e.g, [19]). In general, the VHE emission 
in these objects is explained by the acceleration of charged 
particles in shocks created by the collision of the wind of 
the stellar companion with the expanding pulsar wind or 
by acceleration in accretion-powered jets. The high-energy 
particles produce VHE photons in hadronic and/or leptonic 



(inverse Compton scattering of low-energy stellar photons) 
interactions. 

We present in the following observations of two binary sys- 
tems with VERITAS at energies above 100 GeV: LS I +61 
303 was observed for more than 60 hours during several 
orbital phases between 2008 and 2010 and the Be/X-ray bi- 
nary 1A0535+262 was observed during a giant X-ray out- 
burst in December 2009 for a full orbit. For details on the 
X-ray analysis and results, see [4, 5]. 

2 VERITAS 

VERITAS is an array of four imaging atmospheric- 
Cherenkov telescopes (lACT) located at the Fred Lawrence 
Whipple Observatory in southern Arizona. It combines a 
large effective area ( > 10^ m^) over a wide energy range 
(100 GeV to 30 TeV) with good energy (15-20%) and an- 
gular (^ 0.1°) resolution. The field of view of the VERI- 
TAS telescopes is 3.5°. The high sensitivity of VERITAS 
enables the detection of sources with a flux of 1% of the 
Crab Nebula in about 25 hours of observations. For more 
details on the VERITAS instrument, see e.g. [16]. 

3 LS I +61 303 

The high-mass X-ray binary LS I +61 303 consist of a BO- 
Ve star surrounded by a dense circumstellar disk and a com- 
pact object of unknown nature (neutron star or black hole). 
Optical observations show that the compact object orbits 
the star every 26.5 days [14, 10] on a close orbit, charac- 
terized by a semi-major axis of a few stellar radii only. Pe- 
riodic variability correlated with the orbital cycle has been 
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Figure 1: Average fluxes or upper flux limits per orbital 
phase bin for gamma rays with energies above 300 GeV 
(black points and arrows) and X-rays (blue crosses: Swift; 
blue open circles: RXTE-PCA) from the direction of LS I 
+61 303 as a function of orbital phase observed from Octo- 
ber 2008 to February 2009. The bottom panel shows the re- 
sults averaged over the whole data set; the top panels show 
the results for individual orbits. Upper flux limits (99% 
probability after [15]) are shown for data points with sig- 
nificances less than 3cr (significance calculation after eq. 
17 from [17]). The shaded regions indicate observations 
which had directly overlapping X-ray measurements by ei- 
ther RXTE or Swift. 
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Figure 2: The same as Figure 1, but for the 2009-2010 ob- 
serving season. 
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Figure 3: The same as Figure 1, but for the 2010-2011 ob- 
serving season. 

observed across the electromagnetic spectrum. Long-term 
X-ray observations show that there is additionally variabil- 
ity of the orbital profiles on multiple time scales, from in- 
dividual orbits up to years [23]. Beside the periodic vari- 
ability, intense X-ray flares with flux increases by several 
factors on the time scale of seconds have been observed 
[21]. 

Fermi-LAT observations show that the gamma-ray emis- 
sion at energies about 30 MeV is modulated with the or- 
bital period (peaking after periastron, <l> = 0.225) [1]. 
The energy spectrum shows a sharp exponential cutoff at 6 
GeV. The LAT results reveal overall flux variations not con- 
nected to the orbital movement by up to 40%. In the VHE 
range, LS I +61 303 has been detected by ground-based 
gamma-ray observatories with emission peaking around 
apastron [9, 2]. The measured energy spectra for the apas- 
tron passage can be described by a power law, with a spec- 
tral index of 2.4 and a flux above 300 GeV corresponding 
typically to about 10-15% of the flux from the Crab Nebula. 

Figures 1-3 show fluxes and flux-upper limits per orbital 
phase bin calculated from 64.5 hours of VERITAS obser- 
vations between October 2008 and January 201 1. The cov- 
ered phases by VERITAS observation are influenced by the 
moon cycle which does not allow a complete coverage of 
a single orbital phase^. The strategy for the observing sea- 
son 2008/2009 was to achieve the widest possible coverage 
during an orbital phase with a flux sensitivity in the range 
of 5% of the flux of the Crab Nebula. For 2009/2010, the 
goal was to achieve a deep exposure around apastron in co- 



1. Note that observations taken under moderate moonlight con- 
ditions are included in the presented data set. 
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ordination with Swift XRT observations. In the observing 
season 2010/2011, the phases between superior conjunc- 
tion and periastron were in the focus of the observations. 

These results are remarkable in the following ways: 

• No strong evidence for emission from LS I +61 303 
during the apastron phases was observed during the 
observations taken in 2008-2010. The 2009/2010 
measurements, with a deep exposure of about 17 
hours in the phase interval $=0.5-0.8, result in up- 
per flux limits (99% confidence level) of less than 
5% of the flux of the Crab Nebula. Note that the ob- 
ject has been detected before at flux levels of 10-15% 
of the flux from the Crab Nebula during these orbital 
phases. 

• LS I +61 303 has been detected in September- 
November 2010 close to the phase of superior con- 
junction (much closer to periastron passage), a re- 
gion of the binary orbit previously undetected by 
VHE instruments. A total of 13.9 hours of LS I +61 
303 were accumulated resulting in a post-trials sig- 
nificance of 5.6 standard deviation. The source flux 
above 300 GeV in the $=0-0.1 phase bin was mea- 
sured to about 5% of the flux of the Crab Nebula. 

• The detection of the binary close to superior conjunc- 
tion does not directly implicate that the processes re- 
sponsible for gamma-ray emission above 300 GeV 
are different than those at lower gamma-ray energies. 
The energy spectrum as published by the Fermi LAT 
collaboration with a cutoff at 6 GeV is an average 
spectrum for many orbital periods, the data presented 
here are for a single phase bin at one orbital period. 

The light curves from VERITAS show that beside the vari- 
ability correlated with the orbital cycle, there is additional 
orbit-to-orbit variability with no obvious pattern. The num- 
ber of upper flux limits and gaps in the light curves shown 
in Figures 1-3 indicate the two main difficulties for observ- 
ing LS I +61 303 with current instruments despite the long 
exposure of more than 60 hours: firstly, the orbital period is 
close to the lunar periodicity (lACTs cannot observe during 
periods of bright moonlight) and secondly the fluxes away 
from the apastron regions seem to be close to the sensitiv- 
ity limits of the instruments. Only continuing monitoring 
in the X-ray to VHE energy ranges together with well co- 
ordinated multiwavelength campaigns will allow us in the 
coming years to probe short and long term variability and 
the emission mechanisms more precisely. 
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Figure 4: Swift/BAT counting rate vs time in the 15-50 keV 
energy range. The horizontal line at 0.1 cts/cm^/s indicates 
the trigger threshold for observations with VERITAS. 




Figure 5: Relative orbit of the neutron star around the 
Be star in 1A0535+262. The primary star lies in the fo- 
cus of the ellipse (0,0) and the axis units are multiples of 
the semi-major axis of the orbit. Note that the inclina- 
tion of the system is unknown. Indicated in colors is the 
Swift/BAT counting rate in the 15-50 keV energy range for 
the orbit starting in November 2009. The dashed lines indi- 
cate nights with VERITAS observations, covering the flare, 
apastron, and periods close to periastron. Orbital parame- 
ters after [13]. 



1 A0535+26 is a nearby (distance of ^ 2 kpc) Be/X-ray bi- 
nary consisting of a pulsar of period 103 s in an eccentric 
orbit (e = 0.47) with an 09.7-BO Ille star. The orbital pe- 
riod is 1 1 1 days. Giant outbursts at X-ray wavelengths are 
observed approximately every five years from the direction 



of 1A0535+262. They are thought to be associated with the 
formation of a transient accretion disk [22] . 

VHE gamma rays are thought to be produced through 
hadronic or leptonic scenarios. In hadronic models, protons 
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are accelerated in the magnetosphere of the neutron star 
and impact the transient accretion disk. The pp-interactions 
in the accretion disk produces VHE gamma rays via the 
decay of neutral pions [12]. In leptonic models, electrons 
are acceleration in strong shocks formed by the interaction 
between pulsar wind and stellar wind. These relativistic 
electrons may produce VHE gamma rays by upscattering 
ambient photons. The high-energy photons are expected 
to interact with ambient photons to produce e+e~ pairs, 
leading to orbital modulation of the VHE gamma-ray flux 
[11]. The system is therefore opaque to gamma rays near 
the peak of the outburst, the highest gamma-ray flux is ex- 
pected at the beginning or end of the outburst, when ther- 
mal emission from the accretion disk is at a minimum [20]. 

VERITAS observed 1 A0535-F262 shortly after the onset of 
the giant outburst in December 2009 for a complete orbital 
period. The Swift BAT measurement in Fig 4 shows that 
the hard X-ray flux reached a level of >5 times that of 
the Crab Nebula. The triggering condition for observation 
of flaring X-ray binaries with VERITAS was fulfilled on 
2009 December 5, observations started on 2009 December 
6, shortly after the beginning of the giant flare. They were 
delayed by one day due to very bright moonlight condi- 
tions. The observation covered most of the flare, included 
the apastron phase, and continued for almost 90 days until 
the following periastron phase (see Figure 5). 

No evidence for VHE gamma rays has been found for the 
complete VERITAS data set consisting of 23 hours of high- 
quality observations. The flux upper limit at the 99% con- 
fidence level [15], assuming a power-law-like source spec- 
trum with a spectral index of F = 2.5, is F (> 0.3 TeV) 
< 0.5 X 10-^2 ph cm-2 s-^ (about 0.4% of the flux of the 
Crab Nebula above 0.3 TeV). 

The VERITAS data were arranged in different periods, as 
gamma-ray production and absorption is expected to vary 
with orbital movement and flaring state. The four periods 
are: rising portion and falling portion of the giant flare, 
apastron and periastron. No VHE emission has been de- 
tected in any of these periods; upper limits between 0.9 
and 2.0% of the flux of the Crab Nebula (> 0.3 TeV) have 
been derived. 

X-ray observations with the Swift XRT and RXTE PCA in- 
struments show that the spectra are best fitted with a model 
that consists of blackbody and Comptonized emission from 
thermal electrons at temperatures of approximately 2 keV 
and 6 keV, respectively [4]. 

The non-detection at VHE wavelengths, the results from 
the spectral analysis of the X-ray observations, and the non- 
detection of 1A0535+262 at radio wavelengths indicates 
that there is no significant non-thermal electron population 
in this system. The VERITAS upper limit are placing ad- 
ditionally constrains on hadronic models, with the upper 
limits below the predicted flux levels by e.g. [18]. 

1A0535+262 represents the archetype of the class of Be 
binary systems which exhibit giant outbursts and the VER- 
ITAS observations were the best we could hope to get in 



terms of coverage and exposure. A detection or signifi- 
cantly better flux upper limits will probably have to wait for 
a next-generation ground-based gamma-ray observatory. 

In summary, extensive observations of the two binaries LS 
I +6\ 303 and 1A0535-F262 with VERITAS showed that 
VHE emission processes in binaries are not well under- 
stood. A significant improvement in coverage of the differ- 
ent orbital periods and in sensitivity are probably needed to 
constrain or refute the large variety of theoretical models. 
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